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We perform a preliminary study of the ability of the Higgs decay to four leptons to shed light on 
the top quark Yukawa couplings. In particular we examine whether the h ^ 4i ‘golden channel’ is 
sensitive to the CP properties of the top quark couplings to the Higgs boson. We show that kinematic 
distributions are sensitive to interference of the next-to-leading order electroweak corrections with 
the tree level ZZ contribution. This translates into a sensitivity to the top quark Yukawa couplings 
such that meaningful constraints on their CP properties can begin to be obtained once ~ 300 fb“^ 
of data has been collected at ~ 14 TeV, with significant improvements at higher luminosity or with 
a higher energy hadron collider. This makes the h ^ 4£ channel a useful probe of the top quark 
Yukawa couplings that is qualitatively different from already established searches in h ^ Vy two 
body decays, tth, and gg —>■ h. We also briefly discuss other potential possibilities for probing the 


top Yukawa CP properties in h —>■ 2(.y and t 

INTRODUCTION 

The observation of a Higgs-like resonance with mass 
near 125 GeV [HE] completes the Standard Model (SM) 
and opens up a vast new research program in studying 
its detailed properties in order to determine whether it 
is in fact the SM Higgs. Direct study of the boson itself 
is the best way to unravel the nature of this new state 
and answer interesting questions such as whether its in¬ 
teractions violate CP. It has been established that its 
couplings to ZZ are dominantly CP even [21ll|, but CP 
is violated in nature, so if there is physics beyond the SM 
(BSM), some Higgs couplings may not conserve CP. 

In the SM, the largest coupling of the Higgs is to the 
top quark. Therefore, studying the Higgs top system is 
particularly interesting because it could be an ideal place 
to discover new physics. Furthermore, because of the 
size of this coupling, the hierarchy problem is sharpest 
in the top sector, so potential solutions to the hierar¬ 
chy problem could easily modify the couplings between 
the Higgs and the top. This coupling can be studied di¬ 
rectly using the tth production rate, which as yet is unob¬ 
served |5H7|. Various studies have also shown that kine¬ 
matic observables can be constructed to study the size 
and CP properties of the top Yukawa in this channel 
at the LHC (HHEj) though it requires measurements of 
top and Higgs decays which may be difficult in the high 
luminosity environment of the LHC. 

The Higgs decay to photons is mediated by a top quark 
loop (and the larger W loop), so this channel can also 
probe the top Yukawa coupling. Similarly, the cross sec¬ 
tion of Higgs production via the gluon fusion process is 
sensitive to the top Yukawa coupling. Therefore global 
fits using rates can be used to constrain it under various 
assumptions. This has been done by the experimental 
collaborations [Hdl], as well as by several theoretical 
groups [161425] which indicate a top Yukawa coupling con¬ 
sistent with the SM. Even with the various assumptions 


—>■ hZ, hy. 

in these analyses, sizable deviations from the SM predic¬ 
tion are still allowed. 

Additional probes of the CP properties of the top 
Higgs system include production of a Higgs in association 
with a single top UniESHS], kinematic distributions in 
gluon fusion Higgs production [32l [33] , and low energy 
CP violating observable such as EDMs [33]. All of these 
probes, however, are indirect and suffer from signihcant 
inverse problems. Namely, if deviations from SM predic¬ 
tions are discovered, it is very difficult to determine if 
they are coming from modifications to the top Yukawa 
coupling, or some other type of new physics. Therefore, 
it is important to have as many complementary probes 
as possible. 

In this work we propose a new avenue to study the 
top-Higgs system: the Higgs decay to four leptons. This 
so-called ‘golden channel’ has already been used exten¬ 
sively to study the spin of the Higgs as well as the CP and 
tensor structure of its coupling to gauge boson pairs [3S1- 
EO]. The leading contribution to the golden channel comes 
from the tree level coupling of the Higgs to ZZ generated 
during electroweak symmetry breaking (EWSB). At one 
loop, however, additional couplings of the Higgs to Zy 
and 77 pairs (as well as ZZ) can be generated and me¬ 
diate Higgs decays to four leptons. 

In the SM, these next-to-leading order (NLO) contri¬ 
butions are dominated by W and top loops. While these 
are one-loop contributions, the large available phase 
space for the Zy and yy intermediate states as well as 
the differential spectra allow for these one loop contri¬ 
butions to be distinguished from the tree level ZZ cou¬ 
pling [531 IMl IMI ESI EH]- In particular, due to interfer¬ 
ence effects between the higher dimensional Zy and yy 
couplings with the tree level ZZ coupling, the h ^ 41 
[41 = 2e2/i, 4e, 4/r) channel is surprisingly sensitive to the 
CP properties of these loop induced couplings, especially 
for 77 [661 Eg. 

In this work, we exploit the fact that the top quark 
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mediates Higgs decays to both 77 and Z 7 intermediate 
states via the same couplings to the Higgs boson. There¬ 
fore, unlike previous work which focused on measuring 
higher dimension effective couplings to gauge bosons, we 
here use the underlying loop processes to gain sensitivity 
to the physical parameters of the SM or BSM effects. Fur¬ 
thermore, because the one-loop top mediated effects in¬ 
terfere with the tree-level diagram, the differential cross 
section has a component which is linearly sensitive to 
CP violation and only contains one power of the loop 
factor. This is in contrast to > 77 and h —>■ Zj two 
body decays or tth rate measurements which are sensi¬ 
tive only to the sum of squares of the CP even and odd 
components of the top Yukawa coupling. Therefore, if 
the coupling of the Higgs to the top has the wrong sign, 
as can happen in certain two Higgs doublet [niE 2 ] and 
triplet m models, or if there is a non-trivial CP phase, 
this can in principle be observed directly using only the 
golden channel independently of these other measure¬ 
ments. 

Here we perform an initial feasibility study to explore 
whether the golden channel can be used as a probe of 
the Higgs top quark Yukawa coupling and perhaps un¬ 
cover CP violation. To do this we utilize the parame¬ 
ter extraction framework developed in [Ml IMHMl Ffi] 

to study effective Higgs couplings and adapt it to in¬ 
clude the leading contributions from top quark (and W) 
loop effects. We demonstrate a proof of principle that the 
h ^ M channel has the potential to probe the CP prop¬ 
erties of the top Yukawa at the LHC with very promising 
prospects at a future higher energy hadron collider. We 
also briefly discuss other potential possibilities for prob¬ 
ing the top Yukawa in ft, —>■ 2^7 and hZ, hj. 

PROBING THE TOP YUKAWA IN ft ^ 4^ 

Many previous studies of the golden channel have fo¬ 
cused on probing effective couplings of the Higgs to gauge 
bosons of the form, 

Cvv ~ J Z^^Z, + 

+ , ( 1 ) 

where V,V' = Z,j, and is the usual field 

strength (dual field strength) tensor. These mediate 
Higgs decays to four leptons via the diagram shown 
in Fig. [2 The differential distributions for the many kine¬ 
matic observables in ft —>■ 4£ [sa IMIIMIISZ] give us 
a probe into detailed properties of these effective cou¬ 
plings. In particular, it was demonstrated in [ 66 l [ 68 ] 
that golden channel measurements are surprisingly sen¬ 
sitive to the effective couplings of the Higgs boson to Zj 
and 77 pairs. Because of shape differences and interfer¬ 
ence with the tree level ZZ coupling, the sensitivity is 



the ft — >■ If amplitude where Fi,2 = Z, 7 and £ = e, fi. 

strong enough that SM values of the 77 effective cou¬ 
plings should be probed well before the end of LHC run¬ 
ning. Prospects for Zj are less promising, but still per¬ 
haps possible at the LHC and very promising at a future 
100 TeV collider. This motivates the question of whether 
the sensitivity to these effective couplings translates into 
sensitivity to the underlying loop processes. 

In the SM, the A^^ couplings are generated at one 
loop dominantly through a W boson loop followed by 
the smaller top loop contribution shown in Fig. while 
A^^ is zero at this order. The coupling is generated 
at one loop, but vanishes for an on-shell photon and to 
leading order in the heavy loop particle expansion. The 
leading VF-loop contribution to A^^ involves parameters 
such as the W mass and gauge couplings that are well 
measured from LEP [751175] and the LHC [TJ [T] [771 [TS] 
experiments. Therefore, it is a reasonable approximation 
to take these W loops to be fixed during our parame¬ 
ter extraction of the top Yukawa coupling. Studying the 
sensitivity in ft —>■ 4£ to electroweak parameters in the 
W loops would also be interesting, but requires a more 
careful treatment of other SM one-loop contributions and 
so is left to ongoing work [75] . 

The top loop on the other hand involves various pa¬ 
rameters which are not as precisely constrained. Further¬ 
more, axial couplings between the Higgs and the top 
quark can generate A^"’’ effective couplings which are 
vanishingly small in the SM. Focusing on the top Yukawa, 
we take the top pole mass and the ttZ coupling to be 
fixed, though it would be interesting to study these as 
well [ 8 OI HI] . We parametrize the top Yukawa couplings 
as, 

A 3 + iya^)t, (2) 

where mt is defined to be the pole mass found in the top 
quark propagator with yt = l,yt = 0 at tree level in the 
SM. Note that the large pole mass of the top leads to lit¬ 
tle sensitivity to the top mass in practice in ft —>• VV de¬ 
cays. This is equivalent to saying the top contribution is 
well approximated by a constant effective hVV coupling 
after the top has been integrated out. Thus whether we 
fix or allow the top pole mass to vary makes a negligible 
difference on our results. 
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FIG. 2. One-loop contributions from top quark (left) and W boson to /i —>■ V 1 V 2 —t 4£ (Vi = Z,j). 


After the W and top, the next largest contribution 
to the effective Zj and 77 couplings comes from the 
bottom quark contribution. This effect is suppressed 
by ^ (mb/rrit)^ in the matrix element relative to the 
top contribution which is itself subdominant to the W 
loop. Thus, to a very good approximation, the Z-j and 
77 effective couplings only receive contributions at one- 
loop from the W boson and top quark. 

The h 4£ process receives additional one-loop elec- 
troweak (EW) corrections that are not of the form 
shown in Fig. Since the Zj and 77 effective couplings 
in Eq. ([^ are only first generated at one loop, they do 
not receive a contribution from these additional EW cor¬ 
rections at this loop order. These include processes such 
as corrections to the Z propagator and coupling to lep¬ 
tons as well as various other non-local interactions all of 
which are computable |82[ 183) . Thus in principle we can 
make a precise prediction for all contributions not in¬ 
volving the top Yukawa coupling. This allows us to treat 
this part of the amplitude which does not depend on the 
top Yukawa as part of the SM ‘background’ to our top 
Yukawa ‘signal’. 


Discussion of Signal and ‘Backgrounds’ 


To be more explicit, we can write the > 4£ amplitude 
up to one loop as follows, 

A44£ = ■ (3) 

The leading term arises from the tree level hZZ 

coupling, 


C%M A "^hZf^Z, 


AD 


(4) 


which is generated during EWSB and is responsible for 
giving the Z boson its mass. The second term in¬ 

volves all SM one-loop contributions independent of the 
top Yukawa, though there are one-loop corrections from 
top quark loops to the Z boson propagator for exam¬ 
ple. Finally, Ml encodes the one-loop contribution sensi¬ 
tive to the top Yukawa coupling and which enters via the 
first diagram in Fig. ^ In this work, we will treat Ml as 


our signal and fit for the parameters in Eq. ([^, while we 
will treat the rest of the matrix element as ‘background’ 
which we keep fixed. There are also real non-Higgs back¬ 
grounds, whose leading contributions must be accounted 
for as well and will be discussed below. 

We can further characterize the ‘background’ in M^;^, 
by isolating those contributions which are generated by 
hVV (where VV = ZZ, Z 7 , 77 ) effective couplings of the 
form shown in Fig. [^to write, 

■^EW — -^EW~^-^EW> ( 5 ) 

where we have defined, 

—-^Ew EW ■^'ew (®) 

These contributions all have the form of Fig. and will 
be examined more closely below. 

There are many contributions to M\;yy, all of which 
are computable and can in principle be extracted 
from [HU [53] . Some of these one loop contributions can 
be absorbed into shifts of the tree level couplings. Others 
can be modeled using effective operators. There are also 
real photon emission effects in /i —>■ 4£ [531 - IM] which can 
be non-negligible in certain regions of phase space, but 
which can also be included [85] . The key point however is 
that these corrections do not depend on the top Yukawa, 
allowing us to treat them as fixed when fitting for the top 
Yukawa. Furthermore, since at one loop these corrections 
do not contribute to the Z^ or 77 effective couplings to 
which we are most sensitive in > 4£ |661168) , and since 
they are sub-dominant over most of the phase space [5S] , 
we will neglect them in this preliminary study. However, 
a detailed investigation of their effects is worthwhile and 
will be done in future work. Thus in the end, for the 
present study we define the Higgs part of our ‘back¬ 
ground’ (in contrast to non-Higgs background to be dis¬ 
cussed) as, 

■^BG =-^SM +-^'ew- C^) 

This part of the h ^ 4£ amplitude will be treated as fixed 
during the parameter extraction procedure. 

As mentioned, our ‘signal’ is then the top quark loop 
in the Z'y and 77 effective couplings which we call Mf^ 


1 There is also a wave function renormalization for the Higgs that shapes at one loop and since we are not using the overall rate in 

depends on the top Yukawa, but this does not affect kinematic our likelihood analysis, we can ignore it. 
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and . Of course a top quark loop will also generate 
A4f via hZZ effective couplings, though in practice the 
sensitivity to this contribution is an order of magnitude 
weaker than for the Z 7 and 77 effective couplings | 66 l 
I 68 j . Thus our final signal involving a top loop can be 
written as, 

Ml=M^^+ +Mr- ( 8 ) 

All of the contributions in Eq. 0 and Eq. ([^ enter the 
h ^ A£ amplitude via the hVV couplings and can be 
represented by the diagram in Eig. Thus, by focus¬ 
ing on these and neglecting we are assuming in 

the present study that any deviations from the tree level 
SM prediction occur only through loops which generate 
the hVV effective couplings. As discussed above this is a 
reasonable approximation for current purposes. We will 
examine the contributions in Eq. Q and Eq. ([^ more 
closely below. 

There is of course a non-Higgs background which 
comes dominantly from the continuum qq —>■ M pro¬ 
cess [H] and can have important effects. As discussed 
in [ 68 ] this background enters almost entirely due to 
detector resolution effects. If detectors had perfect en¬ 
ergy resolution, the signal region would essentially be 
a i5-function centered at the Higgs mass leading to an 
effectively background free sample. However, imperfect 
detector resolution has the effect of widening the signal 
region, thus introducing more non-Higgs background into 
the sample and degrading the sensitivity to the hVV ef¬ 
fective couplings [5H] . 

Eor this qq —>■ 4^ background we utilize the analytic 
expressions computed in [56l|64] and follow the procedure 
in gTlISHI to build a signal plus background likelihood 
which includes the parton distribution functions as well 
as crude modeling of detector resolution effects. More 
details on this implementation can be found in [4711561 
IMIEH]- For a more realistic analysis, careful treatment 
of detector resolution and additional background effects 
can be done with the framework in (STHHIESj, but is 
left to future work. However, these detector effects are 
not expected to qualitatively change the results obtained 
here. 


The top and W Loops 

Restricting our attention to the loops which generate 
the hVV effective couplings in Eq. Q and Eq. Q, the 
matrix element for the h ^ Al decay can be written as, 

M{h 4£) = V 1 V 2 ) X 

^ 2£)V,is{V2)M^{V2 ^ 2£) , (9) 

where V = Z, 7 , and 7^^,y(Vi) are the propagators of the 
vector bosons. The second line is described by well mea¬ 
sured physics of vector bosons coupling to leptons, while 


the matrix element on the first line encodes all the Higgs 
physics and for which constraints are far weaker. 

The /i —>■ V 1 V 2 matrix element can be decomposed into 
the following tensor structure, 

^ V 1 V 2 ) = + (10) 

^Clik'(k^ - fci • k2gn + 

where i = .ZZ, ^ 7,77 and ki and k 2 represent the four 
momenta of the intermediate vector bosons (or lepton 
pairs). The Lorentz invariant form factors are in gen¬ 
eral momentum dependent for the off-shell intermediate 
vector bosons and have the generic form, 

Cn--gxMml/mj^,kl,kl), ( 11 ) 


where fi{m\/m\,k\,k\) is the loop function for under¬ 
lying particle X with coupling to the Higgs gx- For a 
~ 125 GeV Higgs mass, the dependance on kf is rather 
weak jSS] over much of the phase space and, to a suf¬ 
ficiently good approximation, the are given by set¬ 
ting kf equal to the physical mass of the relevant gauge 
boson. The kf dependence of the form factors can be 
relevant in certain regions of phase space and factoring 
it in may aid in sensitivity, warranting closer examina¬ 
tion. However, in this initial study we seek to first estab¬ 
lish a proof of principle with the leading terms leaving 
a more detailed exploration of these ‘off-shell’ effects to 
currently ongoing work [75]. 


Thus, the form factors C^g in Eq. (10) will be the ones 
that control Higgs decay to on-shell 77 and Z7 pairs. The 
leading contributions to these form factors comes from 
W and top loops which are shown in Eig [^ These 
one-loop contributions have been computed for h —>■ 
Z7 [87] EH] and h — >■ 77 [89l [90] (including pseudoscalar 
couplings m for the top) and can be straightforwardly 
incorporated into the analytic expressions for the h ^ A£ 
fully differential cross section computed in [35] [M] . For 
our explicit expressions of the top and W loop functions, 
we use the conventions in [92] [93] . 

As discussed above, the sensitivity to the higher di¬ 
mensional hZZ effective couplings is significantly weaker 
than for the hZj and hjj effective couplings [55] ■ Fur¬ 
thermore, though the hZZ effective couplings receive 
contributions from top and W loops, there are also a 
number of other one-loop contributions involving Z and 
Higgs bosons. The already weak sensitivity to these hZZ 
couplings makes disentangling the top contribution from 
other contributions difficult. We therefore simply will 
model these with the set of dimension 5 operators: 


k:zz D + Af^Z'^'Z^, 

+ 4Af^5^Z,Z^") , (12) 


where the are taken as real and constant. To study 
the potential effects of these contributions we treat 
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as nuisance parameters in our parameter extraction pro¬ 
cedure allowing them to vary along with the top quark 
Yukawa. As we will see, the effects of the operators 
in Eq. (12) do not greatly affect our sensitivity to the 
top Yukawa, especially once sufficient statistics are accu¬ 
mulated. 


Other Possible Probes of the Top Yukawa 


distribution with a gaussian of u = 2 GeV centered at 
the Higgs mass which we take to be 125 GeV. Note that 
these resolution effects also enter into the lepton pair in¬ 
variant masses {Mu)- Following the procedure in |47j . 
the parton level differential cross sections for h —>■ 
and qq Ai are combined with the (CTEQ611 [23 HE]) 
parton distributions for the gg and qq initial states. Fur¬ 
ther details and validation of this procedure with Mad- 
Graph5_aMG@NLO [HS] can be found in [MlIBS] . 


In |M| it was shown that due to weak phase/strong 
phase interference effects, the three body h —>• 2i^ decay 
is also sensitive to the CP violation in the effective hZ^ 
and /177 couplings. Thus probing the CP properties of 
the top Yukawa may also be possible in this channel at 
the LHG or future hadron collider. Since this channel is 
less sensitive and requires an understanding of the much 
larger backgrounds than in h —> 41', we do not examine 
this possibility in detail here. 

Grossing symmetry implies —>■ hZ, h'y scattering 

at a future lepton collider jM] |BS] may also be capa¬ 
ble of probing the top Yukawa CP properties. Recently 
it has also been shown that interference between signal 
and background can be used to probe the effective hZ^ 
and /177 couplings in gg —>■ 2^7 |96j . which implies this 
may also be used to probe the top Yukawa. We leave an 
investigation of these interesting possibilities to future 
work. 


SENSITIVITY AT LHC AND BEYOND 

We now quantitatively explore the feasibility of the 
LHG or a future hadron collider to probe the CP prop¬ 
erties of the top Yukawa coupling in h —>■ 41. In partic¬ 
ular, we estimate approximately how many events will 
be needed in h —>• 41 to begin probing values of Yukawa 
couplings which are of the same order as the 0{1) SM 
prediction. We also examine approximately at what point 
h —)■ 41 will become relevant as a measurement relative 
to h —)• V 7 and tth searches for studying the top Yukawa 
(we will not consider gg —> h, but see [161125] for various 
studies of this channel). Once this level of sensitivity is 
reached, a more complete analysis including the various 
other one-loop corrections discussed above will need to 
be conducted in order to give precise constraints on the 
top Yukawa. 

For all results in the present study we have utilized 
the Higgs effective couplings extraction framework devel¬ 
oped in [5B1 IMl IBBl [B7] which incorporates all observables 
available in the (normalized) h ^ Ai fully differential de¬ 
cay width and adapted it to include the top and W loop 
functions discussed above. Also as discussed, we include 
the dominant qq —>■ Ai background and a crude modeling 
of detector resolution |B5j. For the Higgs signal, this in¬ 
cludes a smearing of the four lepton invariant mass (M 4 ^) 


Parameter and Phase Space Definition 


Before presenting our results, we first define our pa¬ 
rameter and phase space. As discussed above, in order to 
study the effects of some of the one-loop contributions we 
have not computed which enter through the ZZ sector, 
we allow the higher dimensional effective ZZ couplings 
in Eq. (12) to vary in the fitting procedure. Thus we de¬ 
fine our multi-dimensional parameter space as. 


^ — (j/t; 




(13) 


Note in particular that we are taking the tree level hZZ 
coupling as fixed and equal to its SM value in Eq. <13- 

To estimate the sensitivity we obtain what we call an 
‘effective’ cf{\) or average error defined in [BS] as, 

^(A) = y|(|A-Ao|), (14) 

where A is the value of the best fit parameter point 
obtained by maximization of the likelihood with re¬ 
spect to A. Here Aq represents the ‘true’ value with 
which our data sets are generated utilizing a Mad- 
Graph5_aMG@NLO [BB] implementation of the effective 
hVV couplings [SBl [M] . The average error is then found 
by conducting a large number of pseudoexperiments with 
a fixed number of events and obtaining a distribution for 
A which will have some spread centered around the aver¬ 
age value. We then translate the width of this distribution 
into our effective cr(A) which converges to the usual inter¬ 
pretation of cr(A) when the distribution for A is perfectly 
gaussian. We repeat this procedure for a range of number 
of signal events (Ns) to obtain a{X) as a function of Ns. 

Following the strategy proposed in [ 68 ], we will use a 
set of phase space cuts which are optimized for sensitivity 
to the Zj and 77 effective couplings. These cuts were 
shown to greatly improve the sensitivity to the Z"/ and 
77 effective couplings over currently used GMS cuts [HBl 
unoj. They are defined as: 


• 115 GeV < Mm < 135 GeV 


• Pt > (20,10, 5, 5) GeV for lepton px ordering, 

• \r]i\ < 2.4 for the lepton rapidity, 

• Mu > A GeV, M«(OSSF) ^ ( 8 . 8 ,10.8) GeV, 
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where Mu are all six lepton pair invariant masses and we 
explicitly remove events with opposite sign same flavor 
(OSSF) lepton pairs that have Mu in the range 8.8 — 
10.8 GeV in order to avoid contamination from T QCD 
resonances. We refer to these as ‘Relaxed—T’ cuts. 

While these cuts perform significantly better in terms 
of sensitivity to the effective hZ^ and /177 couplings than 
the currently used CMS cuts [ 68 ], they also allow more 
non-Higgs background into the sample. It is therefore 
necessary to include the dominant non-Higgs qq —>■ At 
background discussed above as it can have a significant 
effect on parameter extraction when these cuts are uti¬ 
lized. To do this we combine the background and signal 
into a single likelihood and fit for the background fraction 
during the parameter extraction procedure along with 
the parameters in Eq. (13). The background fractions 
used during event generation can be found in [55] . Many 
more details on the various aspects of the parameter ex¬ 
traction framework including the building of the signal 
plus background likelihood and the fitting procedure can 

be found in j??] [SSj [Ml [SSI ISZ] ■ 

We also comment that for these cuts some of one-loop 
EW corrections we have neglected may become 

relevant. For this reason we also will discuss results uti¬ 
lizing CMS-like cuts [ 68 ] for which these contributions 
are phase space suppressed [85], but this will not quali¬ 
tatively affect the discussion. 


Sensitivity as Function of Luminosity 

In Fig. [^we show sensitivity curves for (red) and 
cr(yt) (blue) as function of the number of signal events 
{Ns) (bottom axis) and luminosity x efficiency (top axis) 
assuming SM production {gg —>• h plus VBF at 14 TeV) 
and branching ratios |101[ I102j . In these fits we have uti¬ 
lized the Relaxed—T cuts discussed above and include 
both signal and the dominant qq —> 4£ background. We 
have combined the 2e2p,, 4e, 4/r channels and fit to a ‘true’ 
point of A = (1,0|0.01,0,0.007) corresponding to the SM 
prediction for the top Yukawa which is indicated by the 
dotted black line. 

We see stronger sensitivity to the axial coupling yt 
than to the vector-like coupling yt. This is because the 
CP even component of the top loop is dominated by 
the W loop, but the CP odd couplings yt does not have 
to compete with an analogous W contribution. We also 
study the effect of floating the effective ZZ couplings 
(solid curves) defined in Eq. ( [T^ , versus holding these 
couplings fixed (dashed curves). The values chosen for 
these ZZ effective couplings are only representative and 
whether we take their true value to be zero or 0 ( 10 “^) 
makes negligible difference since the sensitivity to these 
couplings is weak [SZIISH]- What is important to establish 
is whether allowing them to vary in the fit affects the sen¬ 
sitivity to the top Yukawa. We see clearly in Fig. [^that 


this effect is small as expected from differences in the 
kinematic shapes of the ZZ, Z 7 , and 77 intermediate 
states [67l[M] . 


X e(fb-’) 

10 ^ 10 ^^ 



FIG. 3. Sensitivity curves for cr{yt) (top, red) and cr{yt) (bot¬ 
tom, blue) as function of the number of signal events {Ns) 
(bottom axis) and luminosity x efficiency (top axis) assuming 
SM production {gg —>• h plus VBF at 14 TeV) and branching 
ratios [Rmnu^ . In these fits we have utilized the Relaxed—T 
cuts discussed in the text and included both the h ^ Ai 
(4£ = 2e2fi, 4e, 4/r) signal and the qq —>■ At background. We fit 
to a ‘true’ point of A = (1,0|0.01,0, 0.007) corresponding to 
the SM prediction for the top Yukawa which is indicated by 
the dotted black line. We also demonstrate the effect of float¬ 
ing (solid) the effective ZZ couplings (see Eq. lH) versus 
keeping them fixed (dashed). 

The crucial point to emphasize is that we should be 
able to probe 0(1) values of the top Yukawa coupling 
with ^ 6000 — 10000 events corresponding to ^ 800 — 
1500 fb“^ assuming 100% efficiency. Of course in reality 
the efficiency is significantly less, so more realistically 
~ 2000 — 5000 fb“^ may be needed depending on detector 
performance as well production uncertainties. The lower 
ends of this range should be within reach at the high- 
luminosity LHC, and even better sensitivity would be 
achieved with a future hadron collider at higher energy. 


Probing top Yukawa CP Properties 

The results in Fig. [^indicate that the LHC or a future 
collider may be able to directly probe the CP proper¬ 
ties of the top Yukawa coupling in /i —)■ 4.^. To further 
investigate this we show in Fig. [^ and Fig. [^results from 
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c 

g,{tth) 

y{h 77 ) 

y{h Z-i) 

Current 

2.8 ± 1.0 [5] 

1.14 ±0.25 [Ml 

NA 

300 fb'^ 

1.0 ±0.55 [ISi] 

1.0 ± 0.1 [M] 

1.0 ± 0.6 [M] 

3000 fb-i 

1.0 ±0.18 [133] 

1.0 ±0.05 [M] 

1.0 ± 0.2 [M] 


TABLE I. Values of current constraints and future projections 
on the relative signal strength /ii = ojasM (or BR/BRsm) 
for given luminosities. 


the fit for the Icr allowed region in the yt — yt plane for 
a range of data set sizes. The allowed parameter space 
corresponds to the entire region inside the ellipse. 

In addition to utilizing the Relaxed—T cuts (middle, 
yellow ellipses) as in F ig. we also show results us¬ 
ing CMS-like cuts im 11 ( 11 ^ (large, red ellipses). This 
makes it clear the improved sensitivity obtained when 
the Relaxed—T cuts are used. For comparison and as 
a demonstration of the ideal case, we also show the Icr 
region obtained assuming a pure signal sample (inner, 
turquoise ellipses) using these optimized cuts. This also 
makes clear the effects of the qq —>■ 4£ background. 

Fig. i and Fig. also compare the golden channel 
to other measurements which are sensitive to the top 
Yukawa coupling: the tth cross section, the branching 
ratio of —>■ 77 , and the branching ratio of /i —>■ 
Z 7 . The 1 — cr contours are derived from the relative 
signal strength {y,i = ajasM or BR/BRsm) for each 
measurement given by, 

~ _^0.42y2 (15) 

yi{h ^ 77 ) ~ (1.28 - 0.28 j/t)^ + (0.43 y*)" 
y{h ^ Z-i) ~ (1.06 - 0.06 yt)^ + (0.09 y*)", 

where for y,{tth) we use the cross section at 14 TeV for 
the approximate value in terms of y* and yt |105) and the 
numerical factors in /i —>■ V 7 are obtained by evaluating 
the top and W loops [5U [33] at 125 GeV. The values we 
use for the yt signal strengths are summarized in Table |lj 


Before discussing our results further, we comment that 
from the numerical values in Eq. (151, it is clear that the 
sensitivity to the top Yukawa in /i —> 4£ is driven by the 
77 intermediate states. This implies that a reasonable 
approximation of the sensitivity to y* and y* could have 
simply been obtained from a naive rescaling of the results 
for the sensitivity to the 77 effective operators found 
in 1231 |3H|- However, we emphasize that this rescaling 
ignores potential correlations between the Zj and 77 ef¬ 
fective operators [56l|64l|67|. Furthermore, the parameter 
fitting done in this study is qualitatively different since 
(ignoring ZZ couplings) only two parameters {ytiTjt) are 
floated in contrast to four (Hf"'', Hg when us¬ 

ing effective couplings. For these reasons we have not 
simply done a rescaling of the effective couplings, though 
the end results for the sensitivity to yt and yt are not 
drastically different. 


The current Itr confidence intervals obtained in tth 
(green band on the left) [5| and —)■ 77 (blue band on the 
right) |103j direct searches are shown on the left in Fig.|^ 
where 100 —>■ 4^ events have been assumed. We see 
that at this stage h ^ M is not competitive with tth 
and ft, —)■ 77 searches. For 800 events shown on the right 
we use the projected Itr intervals from tth and ft —>■ 77 
searches assuming 300 fb“^ |1041I105] and a SM-like cen¬ 
tral value. We have also added the Icr projections from 
ft —>■ Zj (thick pink band) |106j searches which start to 
become relevant at this luminosity. We can see at this 
stage that ft —>■ 4£ is also starting to become a useful 
channel to complement tth and ft —>■ Y7 searches for 
studying the top Yukawa. 

In Fig. we show the same results, but for 8000 (left) 
and 20k (right) events corresponding to > 1000 — 3000 
fb“^ and where the projected Icr intervals from tth, 
ft —)■ 77, and ft —)■ Z^ searches have been used assuming 
3000 fb“^ |104H1()6] . We see in these results that if we as¬ 
sume the Higgs couplings to ZZ and WW are positive, 
eventually ft —>■ 4£ should be able to establish the overall 
sign of yt independently of any other measurements of 
the top Yukawa. We further see the possibility of using 
ft —4^ as a consistency check with tth and h ^ Vj 
searches as well as the qualitatively different nature of 
the ft —>■ 4^ measurement. 

The results in Fig. and Fig. make it clear that 
ft —>■ 4 .^ is a useful and complementary channel to tth, 
ft —>■ Z^, and ft — 77 searches for probing the top 
Yukawa at the LHC or a future collider. Furthermore, 
depending on how sensitivities evolve over time, it may 
be possible that h ^ 41 will be able to constrain re¬ 
gions of parameter space which are difficult to probe in 
other channels helping to ensure that potential CP vio¬ 
lating effects would not go unnoticed. In the event where 
a deviation from the SM value is observed in either on- 
shell ft —>■ Zj, 77 two body decays or tth production, 
the four lepton channel will be a crucial ingredient in 
both confirming and characterizing the anomaly. Quan¬ 
tifying more precisely these possibilities will require a 
detailed treatment of the various one-loop and off-shell 
effects which we have not included, but a thorough inves¬ 
tigation is left to ongoing work m- Many more results 
from the current analysis can be found in |107j . 


CONCLUSIONS 

We have demonstrated that the ft —4^ ‘golden chan¬ 
nel’ can be a useful probe of the top Yukawa at the LHC 
and future colliders. We have considered the leading ef¬ 
fects in order to give a proof of principle that this channel 
can serve as a complementary, but qualitatively different, 
measurement to ft —)■ 77 and ft —)■ Zj two body decays 
as well as gg —i' h and tth searches for studying the top 
Yukawa. A detailed study of the sub-dominant one-loop 
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FIG. 4. Left: la contours for yt vs. yt with 100 h ^ 4£ events corresponding to ~ 15 — 40 fb~^ at the LHC14 assuming SM 
production and branding fractions | 10111102| and depending on detector efficiencies. The allowed parameter space is the entire 
region inside the ellipses. The same fit as in Fig. with floating ZZ couplings is performed with the true point represented 
by the star and corresponding roughly to the SM prediction. We show the la confidence interval obtained in h —>■ 4f utilizing 
CMS-like cuts |86II100 | (large, red ellipse) and compare it to the Relaxed—T cuts (middle, yellow ellipse) described in text and 
introduced in m- For comparison with the ideal case we also show the projected la interval assuming a pure signal sample 
(small, turquoise ellipse) and utilizing the Relaxed—T cuts. The current la conhdence intervals obtained in tth (green band 
on the left) [5] and /i —> 77 (blue band on the right) [103] direct searches are also shown (see Table |I|. Right: Same as left, 
but for 800 h —>■ 4f events corresponding to ~ 100 — 300 fb“^. The projected la intervals from tth and h —>■ 77 searches have 
been used assuming 300 fb“^ [10411105] . We have also added the la projections from h —>■ Z 7 (thick pink band) [106] searches 
which start to become relevant at this luminosity. 


and ofF-shell effects in order to quantify the sensitivity to 
the top Yukawa more precisely is ongoing. 

In particular the h —>■ 4f channel can be used to di¬ 
rectly study the CP properties of the top Yukawa in a sin¬ 
gle channel independent of other measurements. This is 
useful because multiple measurements need not be com¬ 
bined allowing us to avoid complications from combining 
errors in different channels in order to establish the CP 
properties. Furthermore, the experimentally clean nature 
and high precision with which this channel is measured 
along with the fact that it is theoretically very well un¬ 
derstood makes it valuable as both a consistency check 
for other channels as well as perhaps the most direct way 
to uncover potential CP violation in the top Yukawa. 

The main drawback of h —)■ 4£ is that it is statistics 
limited, but our results indicate that the necessary preci¬ 
sion to begin probing the top Yukawa may be reached at 
the LHC and certainly at a future hadron collider. The 
theoretical importance of the top Yukawa coupling has 
been firmly established for quite some time and finding as 
many independent probes to study it will be crucial. We 
thus encourage experimentalists to add h —4£ to the 
list of already established channels for studying the top 
Yukawa and in particular its CP properties. 
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